Neutrino Signatures on the High Transmission Regions of the Lyman-a Forest 
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We quantify the influence of massive neutrinos on the statistics of low density regions in the 
intergalactic medium as probed by the Lyman-a forest at redshifts z = 2.2-4. Based on mock 
quasar spectra from hydrodynamic simulations with cold dark matter, baryons and neutrinos, we 
find that the probability distribution of weak Lyman-a absorption features, as sampled by Lyman- 
a flux regions at high transmissivity, is strongly affected by the presence of massive neutrinos, 
when the amplitude of the matter power spectrum is fixed on large scales from measurements of 
the microwave background anisotropies. The application of the derived statistic to observed quasar 
spectra would provide a new sensitive measure of neutrino masses. 



Introduction. Neutrino oscillation experiments re- 
vealed that neutrinos are not massless particles. Since 
then a major effort has been dedicated to measure or 
constrain neutrino masses. Current laboratory bounds 
constrain the electron neutrino mass to m l/c < 2.05 eV 
[l|, 0] • Cosmological bounds for the sum of all neutrino 
masses are still significantly stronger: constraints from 
WMAP7 alone yield Ejm^ < 1.3 eV 0, while combined 
with large scale structure (LSS) measurements they con- 
straint the mass to S^m^ < 0.3 eV 0-0|- The tightest 2a 
upper limit of Ynm Vi < 0.17 eV, is obtained by combining 
cosmic microwave background (CMB) results, LSS and 
Lyman-a forest Q data sets (see [i| for a summary of cur- 
rent and future neutrino mass constrains) . Among all the 
different observables the Lyman-a forest is particularly 
constraining since it probes structures over a wide range 
of redshift, in a mildly non-linear regime and on small 
scales where the neutrino signature is present lol- 12|. 

The dynamics of cosmological neutrinos is very differ- 
ent from that of the dominant cold dark matter (CDM) 
component. The large velocity dispersion of neutrinos 
suppresses their power spectrum of density fluctuations 
on small scales, making the shape of the total power spec- 
trum a potential probe of neutrino masses. 

Previous studies have addressed the role of neutri- 
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the intracluster medium |22j, using both linear theory 



LSS 1642J and 



and N-body /hydrodynamic techniques for the non-linear 
regime. It has been shown that on scales of 1-10 /i -1 Mpc 
the non-linear suppression is redshift and mass depen- 
dent in a way that is different from a naive extrapolation 
of linear theory. 

In this Letter we study the effect of massive neutrinos 
on the properties of low density regions or voids in the 
intergalactic medium (IGM). Naturally, neutrinos have 
only a mild effect on dark matter halos (I34l5( . since 
their large velocity dispersion prevents their clustering 
on small scales. In contrast, we find that the impact of 
neutrinos on void properties is much stronger. Voids are 
relatively empty regions with 6 = p m /p m —l ranging from 
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FIG. 1. Linear power spectra for different neutrino masses 
at z = (upper lines) and z = 3 (bottom lines). The inner 
panel shows the linear (solid lines) and non-linear (dotted 
lines) power spectrum at z=3 normalized to the case without 
neutrinos. On scales k < 0.03 h Mpc^ 1 changes in the non- 
linear power spectra are driven by the differences in the linear 
power spectra. 



almost —1 in their cores to ~ —0.7 at radii 10 — 20 Mpc 
at z = [23j]. By solving the dynamical equations for an 
isolated spherical top-hat undcrdense perturbation, we 
find that neutrinos modify the evolution of undcrdense 
regions by making them smaller and denser. Neutrinos 
contribute to the interior mass of the undcrdense region 
delaying the rate at which CDM is being evacuated from 
its interior and slowing down the velocity of the shell 
surrounding it. We find that the linearly extrapolated 
density contrast when the undcrdense region enters into 
its non-linear phase decreases by ~ 10% for neutrinos 
with Yiim Vi ~ 1 eV. Using the analytic model presented 
in 24| we find that the statistics of voids depend on both 
a 8 and Y>im Vi . Lyman-a voids and their dependence on 
other cosmological parameters have been studied in (2{| . 



Here we study the dependence of void properties on the 
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sum of the neutrino masses. We focus our attention on 
the Lyman-a signature of low density regions, and in- 
troduce a new and simple statistical tool that samples 
most of the IGM volume and appears highly sensitive to 
neutrino masses. 

Numerical Method. Our mock quasar spectra are based 
on cosmological simulations run with the TreePM-SPH 
code GADGET-3 [H]. The code has been extended to 
include neutrinos either by solving their potential on the 
mesh or by representing them as discrete particles 21 



Here, we use primaril y th e first implementation and refer 
the reader to 17|, [l8j, l2l| for a critical comparison of the 
two methods. Our simulations consist of 2 x 512 3 CDM 
plus gas particles sampling a periodic box of 512 /i~ 1 Mpc. 
We adopt a flat ACDM background with cosmological 
parameters Ocdm + fl v = 0.25, Ha = 0.7, S!;, = 0.05, 
h = 0.7 and n s = 1. We consider three degenerates neu- 
trino species with a total mass of E^m^ = 0.0, 0.3 and 0.6 
cV. The initial power spectra of most of our simulations, 
produced with CAMELj, are normalized for all neutrino 
masses at a wavenumbcr 2 x 10~ 3 h Mpc -1 , correspond- 
ing to the scale constrained by CMB data. This produces 
different values of er 8 = 0.877, 0.806 and 0.732 at z = 
for the models with X^m^ = 0.0, 0.3 and 0.6 eV, respec- 
tively (see Fig. [TJ. Our initial conditions are generated 
at z = 49. 

For each simulation we consider snapshots at redshifts 
z = 2.2 and z = 4 that bracket the range of interest 
for the observed Lyman-a forest in quasar spectra from 
ground-based telescopes. For each snapshot we sam- 
ple 4500 random line of sights (RLOSs) uniformly dis- 
tributed along each x, y or z direction. For each RLOS we 
extract the baryon density contrast PhM/Ph an d the pe- 
culiar velocity V p (r) along the line of sight and then com- 
pute the transmitted flux e~ T ^ in rcdshift space (with 
u in km s _1 ), where t is the Lyman-a optical depth, by 
using the Fluctuating Gunn Peterson Approximation: 



r(u) = A [ °° dx S[u -x- VJx)} ( 
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where x = H{z)r/ (1 + z) is the rcdshift space coordinate 
and A is a factor that depends on the global thermal 
history of the IGM [27j], 
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with Tjji being the hydrogen photoionization rate. The 
power-law index in the scaling with Pb/p~b arises from 
the equation of state for the IGM temperature, T = 
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FIG. 2. Real space distribution of the baryon density con- 
trast, Pb/~P~b (t°P panel), and the peculiar velocity, V p (mid- 
dle), along a random line-of-sight (RLOS). In the bottom 
panel we plot the transmitted flux F — e~ T , in units of the 
mean flux (F), in redshift space. 



T (pb/p b ) a [28|, with a ps 0.6. In all our calculations 
we adopt To = 10 4 K and choose Thi such that the 
mean flux over the whole set of RLOS reproduce the ob- 
served mean flux at redshift zJM (F) = e _Te «W with 
T e ff(z) = 0.0023(1 + z) 3 - 65 [30^ We neglect effects of 
thermal broadening. Finally, we smooth the flux over a 
scale of 1 /i -1 Mpc which is larger than the Jeans length 
to avoid sensitivity to substructure below the Jeans scale 
which is affected by numerical resolution and astrophys- 
ical processes (e.g. feedback from galactic winds). 

Figure [5] shows the baryon density contrast, Pb/Pb, and 
peculiar velocity, V p , extracted along a RLOS as a func- 
tion of the comoving coordinate r together with the cor- 
responding transmitted flux F = eT T in redshift space, 
plotted in terms of the mean flux at redshift z. 

Analysis of the simulations. We focus our analysis on 
the statistical properties of low density regions that pro- 
duce weak absorption features. A region is a continuous 
domain in the transmitted flux profile which remains al- 
ways above a given threshold. The higher the threshold, 
the lower the absorption in that region. For each RLOS 
we extract the transmitted flux from Pbl~Pb an d V p and 
count the number of regions above the selected threshold. 
This results in a statistical estimate of the low absorp- 
tion contribution to the Lyman-a signal, and allows us 
to quantify the impact of neutrinos on those regions. 

In Fig. |3] we plot the probability distribution func- 
tion (PDF) for the number of regions per path length 
of 100 ft- 1 Mpc0 above a threshold of F/(F) = 1.14 at 
redshift z = 2.2 (top) and at redshift z = 4.0 for a thresh- 
old F/(F) = 1.70 (bottom) for three different neutrino 
masses, S^m^ = 0.0, 0.3, 0.6 eV. We have verified that 
these PDFs do not change if we increase the number of 
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2 Non-integer numbers are due to the path length normalization. 
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FIG. 3. Probability distribution function (PDF) for the num- 
ber of regions per path length of 100 /i _1 Mpc above a thresh- 
old of F/(F) = 1.14 (top), 1.70 (bottom) as a function of 
Ejiri,/; and erg at z = 2.2 (top) and 2 = 4 (bottom). The 
PDFs have long tails with a very low probability that extend 
up to 10-12. The as — fi„ degeneracy is not perfect and can 
be broken by studying the spectra at different redshifts. 



RLOS, i.e. our statistical sample of RLOS is large enough 
to reliably measure the PDF. Figure [3] shows that the 
neutrino mass has a significant impact on the mean of the 
distributions. In Fig. @] we plot the mean of the distribu- 
tions, i.e. the average number of regions per path length 
of 100 /i _1 Mpc above a given threshold, as a function 
of the threshold for the three different neutrino masses 
(Sim^ = 0.0, 0.3, 0.6 eV) at redshift z = 2.2 (top) and 
z = 4.0 (bottom). This shows clearly that the higher the 
threshold, the larger are the differences between the var- 
ious neutrino cosmologies. This is the expected neutrino 
signature as we discuss below. 

<r%-£l v degeneracy. In order to investigate the degen- 
eracy between <j$ and fi„ we run one simulation with 
SjTO^ = 0.6 eV and the same a% = 0.877 as the 
Him Vi = 0.0 eV model. This mimics the case for which 
erg is measured independently and one would like to con- 
strain fi„. We find that the PDF for the number of re- 
gions above a threshold of 1.14 at z = 2.2 is close to the 
one with Tnm Ui = 0.0 eV (see Fig. [3]), but not identi- 
cal. Larger differences between these two models show 
up at z = 4.0 (see Fig. |4]). Given that neutrinos affect 
the growth factor, we conclude that the degeneracy be- 
tween o"8 and fi„ can be broken by examining the redshift 



evolution of the PDFs. 

Numerical convergence. We explicitly checked that rel- 
ative differences between our neutrino models are nu- 
merically converged against mass and spatial resolution. 
Furthermore, wc used the neutrino particle implemen- 
tation and found the same trends in the neutrino sig- 
nature as with the grid method, although relative differ- 
ences between the different models are even slightly larger 
when we use the particle implementation. This is due to 
the fact that non-linear neutrino effects, such as phase 
mixing, are only captured properly by using the parti- 
cle implementation. We note however that the grid im- 
plementation in the mildly non-linear Lyman-a regime 
is fully justified since non-linear neutrinos effects should 
not be particularly important at those redshifts and at 
k < 1 h Mpc -1 . 

Discussion and Conclusions. We have presented a 
novel and simple method to measure neutrino masses 
by studying the properties of low density regions in the 
Lyman-cv forest of quasar spectra. Those regions corre- 
spond to the inner parts of non-linear voids. We find 
that the number of regions above a given threshold in 
the flux is strongly affected by neutrinos, especially once 
the amplitude of the matter power spectrum is fixed at 
large scales (i.e. normalized by the CMB anisotropies) . 
The changes between different models are due to two 
factors: the change in amplitude and slope in the lin- 
ear power spectrum driven by neutrinos and non-linear 
effects associated with CDM and neutrinos (note that 
neutrinos modify the non-linear evolution of the CDM 
distribution). The inner panel of Fig. [1] shows the linear 
(solid lines) and non-linear (dotted lines) versions of the 
power spectrum at z = 3 normalized to the case with- 
out neutrinos. Whereas the modification on large scales 
(k < 0.03 h Mpc -1 ) is driven by the linear power spec- 
trum, we find that on smaller spatial scales the non-linear 
effects dominate. We have also shown that the as — fit- 
degeneracy can be broken by studying the redshift evo- 
lution of low density regions in the Lyman-Q forest flux. 

Our method can be used to constrain neutrino masses 
by measuring the average number of regions as a function 
of threshold. In the subplot of Fig. 0] we show the av- 
erage number of regions per path length of 100 /i -1 Mpc 
as a function of the threshold at z = 2.2 normalized to 
the neutrinoless model. The black error bars show the 
90% (interior tick marks) and 99% (exterior tick marks) 
confident intervals for a mock catalog consisting of 200 
RLOS taken from the simulation with (E^m^ = 0.0 eV, 
erg = 0.877). We find that with a catalog consisting of 
200 RLOS we can rule out models (£im„ 4 = 0.3 eV, 
(7 8 = 0.806) and (S^m^ = 0.6 eV, cr 8 = 0.732) with 
a high significance. The best way to distinguish mod- 
els with the same as would be by adding more quasar 
spectra and combining results at different redshifts. Al- 
though regions at high transmissivity are prone to sys- 
tematic errors (such as the continuum fitting procedure 
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FIG. 4. Average number of regions per path length of 
100 /i _1 Mpc as a function of flux threshold at redshift z = 2.2 
(top) and 2 = 4 (bottom) for different neutrino masses and erg. 
The subplot in the upper panel shows the ratio between mod- 
els with T,im Vi 7^ 0.0 and the model with Ejm^ = 0.0. The 
black error bars indicate the 90% (interior tick marks) and 
99% (exterior tick marks) confident intervals for a mock cat- 
alog consisting of 200 RLOS taken from the simulation with 
{T,im Vi = 0.0 eV, a 8 = 0.877). Models with E»m^ = 0.3,0.6 
and as = 0.806,0.732 respectively can be ruled out with a 
high significance by using a catalog of 200 QSO spectra. 

which provides the main systematic uncertainty in cal- 
ibrating the flux at weak absorption levels), we believe 
that statistics of large-scale voids such as the one pre- 
sented here could be important for upcoming and present 
spectroscopic surveys of quasars (e.g. [3lj). Variations in 
the thermal and ionization histories of the IGM are ex- 
pected to impact on the Lyman-a properties of large size 
voids j2f| and the neutrino signatures should be seeked 
by marginalizing over all the other relevant parameters. 
However, the generic redshift evolution and threshold 
(scale) dependence of the neutrino effects can be used 
to separate them from the redshift evolution of the IGM 
properties. Independent constraints on thermal evolu- 
tion obtained by different techniques and data sets can 
be put as priors in the analysis performed. 
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